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Fundamental inconsistencies of superstatistics are highlighted. There is no such thing as a su- 
perposition of Boltzmann factors; what is actually derived is a generating function and not a nor- 
malizable probability density. The beta density of the second kind for the energy is shown to be 
subordinated to x^-probability distribution, and this provides a power-tail distribution for large val- 
ues of the energy, comparable to a Pareto distribution. The authors' claim that the Laplace method 
will only work on the inverse Laplace transform in the opposite limit of small energies is shown to 
be incorrect. The exact expression is obtained by contour integration at a multipole located at the 
origin, corresponding to the high temperature limit. 
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Notwithstanding the criticisms logged against 
'superstatistics' Q, the manuscipt 'Asymptotics of 
Superstatistics' has recently been posted on this site 0- 
The basic idea is to obtain 'fat' tail power distributions 
by a 'mixing' of Boltzmann factors. The rationale behind 
why such a mixing of Boltzmann factors will produce fat 
tail power laws is left to the reader's intuition. 

A 'mixing' distribution, is used to obtain a new 

'distribution', B{E), through a Laplace transform 
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where E is the energy and f3 is its conjugate intensive 
parameter. It is claimed that both and B{E) are 
unimodal. This is unnecessary and misleading: all that is 
required is that they be monotone In the manuscript 
under discussion, the Laplace integral, Q was evaluated 
in the high energy, E — > oo, /? — > limit, using Laplace's 
method. Yet, when the inverse transform is considered. 
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along the vertical line = 6, in the region of existence of 
B{E), we are surprised to learn that the same method can 
only be applied in the opposite limit as E 0, P oo. 
This appears all the more disconcerting since no informa- 
tion is lost in a Laplace transformation: A distribution is 
uniquely determined by its Laplace transform and vice- 
versa. Both should yield the same expressions in the 
same asymptotic limit. 

Historically, any relation giving the asymptotic behav- 
ior of B{E) in terms of is known as a Tauberian 
theorem, whereas the inverse relation is referred to as 
Abelian The behavior of B{E) as — > oo deter- 
mines completely the behavior of near the origin 
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and vice-versa. Without loss of generality, we can con- 
sider power laws for both and B{E), one of which 
must be monotonically increasing, while the other mono- 
tonically decreasing. That the Laplace transform ex- 
ist, it is necessary that the exponent in 
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be confined to the semi-open interval [0,oo). This then 
gives 



B{E) = E- 
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which is the tail distribution that Touchette and Beck de- 
rive by Laplace's method in the limit as i? — > oo. They 
are quick to point out that |0J has the same asymptotic 
form as Tsallis' g-exponential distributions 4], where 
m — {q — 

However, the Laplace transform |^ defines a generat- 
ing function, B{E), and not a normalizable probability 
density. The conjugate distribution is the distribution 
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characterized by the 'parameter', i?o, and not the prob- 
ability distribution Q 



where 



p{E) = B{E)/Z, 



B{E) dE, 



since the integral diverges. Hence, a simple Laplace 
transform on the 'density of states', /(/?), will not lead 
to a power-tail distribution. 

The x^-distribution jSJ for the inverse temperature 
can, however, be derived from the x^"distribution for 
the energy by Bayes' theorem of inverse probability j^. 
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Moreover, it is the distribution of the directing process 
which transforms the x^-distribution 
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for the energy into a beta distribution of the second kind 
for the energy. This is accompHshed through subordina- 
tion in which the parameter, /So, referred to as the 
'state of nature', and characterizes the heat reservoir jis 
replaced by the lower hmit on the energy distribution [3 . 

The two x^-distributions, |SJl and are related by 
the Levy transformation 

(3qE = (3Eo, 

where the subscripted quantities refer to the heat reser- 
voir. In analogy with Brownian motion, where we can 
consider a distribution of particles at a given instant in 
time or, equivalently, consider a single particle and fol- 
lows its motion in time, the conjugate variables of posi- 
tion and time become those of energy and temperature 
in thermodynamics The 'pure' spatial process, where 
the initial position replaces the time parameter, is the 
Cauchy distribution, and we say that the Cauchy process 
is subordinated to Brownian motion 0, p. 348]. Analo- 
gously, in thermodynamics, the beta density of the sec- 
ond kind is subordinate to the x^-distribution ©, since 
0, p. 56] 
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where •) is the beta function. For Eq <C E, the beta 
density Q transforms into the Pareto density 
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The energy Eq is then interpreted as the smallest en- 
ergy for which the distribution is valid, corresponding to 
the minimum income in Pareto's distribution of income 
above a given value. 

Touchette and Beck claim that Laplace's method will 
only work on the inverse transform (j^J in the limit as 
/3 — > cxD and i? — > 0, i.e., the low temperature limit. The 
distribution they find by using Laplace's method is 
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where the prime means differentiation with respect to the 
argument which is the implicit solution to 



{lnB{E)y ^~f3. 
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However, PE+lnB{E) has a minimum at E[3 = m/ f3 and 
not a maximum, as contended by the authors. Hence, 



Laplace's method will not work regardless of which 
asymptotic limit is being considered. It is also surprising 
that the authors find it necessary to work in the opposite 
asymptotic limit when discussing the inverse transform. 
No information is lost upon taking the Laplace transform 
or its inverse so the two should necessarily correspond to 
the same asymptotic limit. 

As is well-known, the inverse Laplace transform ^ can 
be performed by contour integration. There is a multi- 
pole of order m at the origin, and B{E) ^ a.s E ^ 00. 
If we enclose the pole at the origin by a contour Cm, we 
have from Cauchy's theorem 
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which is precisely the density ||2J) that we started out 
with. 

Having come this far, we may like to inquire into the 
physical meaning of and B{E). According to the 
process of subordination, the primary process is the 
distribution ® for the energy. The directing process 
is governed by the inverse temperature distribution ((SJ. 
This is analogous to randomizing time in Brownian mo- 
tion where the Brownian motion process for the incre- 
ments of a Brownian particle is transformed into a Levy 
distribution for the random increments in time. 

Equilibrium statistical mechanics teaches us that the 
fundamental quantity is the im^groper probability distri- 
bution, or 'structure' function 
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representing the density of states in phase space. The 
conjugate probability distribution is the exponential law 
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where the generating function is 



n{E), 



Z{(3) 



n{E)e 



-pE 



[3- 



(10) 



(11) 



The functions /(/3), Q and B{E), Q refer to the in- 
verse of the generating function, (|11() , and the inverse of 
the structure function, (jS)), respectively. Hence, they are 
devoid of any physical content. 

Khinchin 8] [p. 87] contends that the conjugate dis- 
tribution (|10|l can be well approximated by the normal 
distribution on the strength of the central limit theorem. 
Setting 
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and evaluating it at (InZ)' 
result 



-/? gives the important 
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where the mean energy, E = m/f] is given by the high 
temperature, equipartition limit. Now, introducing this 
explicitly into ((T^ gives (i. p. 171] 
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which is precisely Stirling's expression for m!. 

The reason for this is that the normal law H13|l applies 
in the limit where the number of degrees of freedom of 
the system has been allowed to increase without limit. 
Rather, if we maintain a finite number of degrees of free- 
dom, and identify the conjugate distribution with the 
distribution ^ for the energy, then we obtain the ex- 
pression for the density of states as 
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without any approximation. Obviously, if we were to 
evaluate (fnji at E = E ^ m/P, we would again get 

/ — +- - 
Stirling's approximation, ml = v27rm'" 2e ™. 

It is appropriate, at this point, to ask for the distinc- 
tion between the x^" ^nd normal distributions as far as 
the processes are involved. This is a matter for order 
statistics. The probability that one of the n — 1 molecules 
has energy falling in the interval from Ei to Ei + dEi is 
(n — \)dEi/ E. There are then n — 2 molecules remaining, 
TO — 1 of which have energies < dEi while n — to — 1 have 



energies > dEi . The probability that to — 1 have energies 
< El is {El/ E)™-~^ , while those having energies > Ei is 
(1 — Ei/Ey^ . The product of the three factors is 
proportional to the beta density of the first kind: 
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The mean energy of the molecules with energies < Ei 
is simply the fraction (m/n) of the total energy. If n is 
so large that we may use the equipartition result for the 
total energy, E = n/ (3, and let n — > oo, the beta density 
of the first kind will transform into the ^^-distribution. 
The mean value Ei ~ m/ P is again given by the 
equipartition result, which is entirely independent of the, 
asymptotically infinite, total number of particles, n. 
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